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^ ! Abstract 

The Hartree-Fock-Bogolyubov (HFB) problem for the cutoff local energy-density 
functional is solved numerically by using the Gor'kov formalism with an exact treat- 
ment of the particle continuum. The contributions from the resonant and "gas" 
continuum to the spectral density of the HFB eigenstates as well as the shifting 
and broadening of the discrete HF hole orbitals are clearly demonstrated with the 
illustrative example of the drip- line nucleus 70 Ca. The structure of the neutron den- 
sity distribution in the localized ground state is analyzed, and the formation of its 
"V ! extended tail ("halo") is shown to be a collective pairing effect. 
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Pairing correlations significantly influence the properties of finite nuclei pro- 
ducing, e.g., the observed odd-even staggering of the binding energies and 
radii. In the mean-field models, the pairing is introduced either at the BCS 
or the HFB level [1]. The zone of the active phase space around the Fermi 
surface, in which the major pairing effects are developed, depends on the ef- 
fective interaction used. Whenever this zone includes states from the particle 
continuum, and this is obviously inevitable when approaching the drip line, 
the BCS approximation breaks down, yielding unphysical nucleonic gas [2] . A 
correct treatment of the pairing problem with continuum is provided by the 
HFB theory which, for the chemical potential /i < 0, always gives a localized 
ground state [3]. In principle, an exact solution of the HFB equations should 



1 E-mail address: fayans@mbslab.kiae.ru 

2 Supported by the Deutsche Forschungsgemeinschaft 



Preprint submitted to Elsevier Preprint 8 February 2008 



be found with the physical boundary conditions both for the bound and scat- 
tering states (examples of such solutions are given in [4,5]). However, in most 
applications of the HFB theory to nuclei, with different kinds of the energy- 
density functionals, the particle continuum is discretized in a box (see, e.g., 
[6,7] and references therein), and it is still questionable whether or not the 
particle level density can be well reproduced with such a prescription, partic- 
ularly near the threshold. Continuum effects have been studied very recently 
in [8] within the HFB method using a fixed, analytically soluble mean field 
potential and keeping self-consistency only in the pairing channel; again, the 
continuum states have been discretized in a spatial box. An attempt to take 
into account the resonant continuum within the HF+BCS approach has been 
also made recently in [9,10]. One should notice that, for a nucleus with even N 
(we consider the neutron subsystem for definiteness) , all HFB solutions with 
energies E > \/i\ rs S2 n /2 (S^n is the two-neutron separation energy) belong to 
the continuum, and these solutions contain both localized and non-localized 
HFB resonances. Their positions and widths, as well as the genetic relation 
to the Hartree-Fock (HF) single-particle spectrum, and also the contribution 
from the nonresonant continuum are worth to be revealed with numerically 
"exact" HFB calculations. Of particular interest is the situation on the drip 
line where the nuclei have only one bound state — the ground state, and the 
whole HFB spectrum is continuous. Our study is motivated by the great ac- 
tivity in physics of radioactive nuclear beams which have lead already, e.g., to 
the discovery of halo structure in some loosely bound nuclei. 

Our approach is based on the generalized variational principle applied to the 
cutoff (local) energy-density functional [11-13] and on the physically trans- 
parent and mathematically elegant Green's function method. The coordinate- 
space technique used to solve the Gor'kov equations exactly, for spherical sys- 
tems with a local pairing field, is described in detail in [4,11,13]. Here we shall 
discuss the results for Ca isotopes obtained with functional DF3 whose normal 
part is parametrized in Ref. [14]. In the anomalous (pairing) part we employ 
the contact "gradient" pairing force [15] with a density-dependent form factor 

J* = 1.35 x C (ft + h^x 2 '\r) + f%rl (Vs(f 



where x = (p n + p p )/2p is the isoscalar dimensionless density with p n ( p ) the 
neutron (proton) density and po the equilibrium density of one kind of particles 
in symmetric nuclear matter with the Fermi energy e F! Co = 2e p/3p is 
the inverse density of states at e = e 0F ; r^ = (3/87rp ) 2//3 (for DF3, C = 
308.2 MeV-fm 3 and r = 1.147 fm). The pairing strength parameters /| x = 
— 1.79, ht = 1.36 and /|? = 2.0 have been deduced from the neutron separation 
energies and charge radii in the Pb isotope chain [13], both the staggering 
and kink being well reproduced. As shown in Ref. [13], the force (1) with 
the scaling factor of 1.35 allows to nicely reproduce the analogous differential 



observables in the Ca isotopes including the anomalous behavior of charge 
radii. All calculations are done with a fixed energy cutoff of 40 MeV above the 
Fermi level in any given nucleus. 

The HFB problem with the above energy-density functional is solved itera- 
tively by using the coordinate-space technique in the complex energy plane. 
The normal and anomalous densities at each iteration are evaluated as the con- 
tour integrals of the generalized Green's function G which obey the Gor'kov 
equation (E — 7i)G(E) = I with 7i the HFB Hamiltonian. For spherical sys- 
tems, after separating the spin-angular variables one gets the equation for the 
radial component cjji of the generalized Green's function G: 
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where h jt = ^(-^ + ^) + U c (r) + U^r)^)^ is the Hartree-Fock (HF) 
Hamiltonian in the jl channel with U c and U s i the central and spin-orbit poten- 
tials, respectively; (al)ji = j(j + 1) — 1(1 + 1) — |; E > 0. The solution of this 
matrix equation is constructed by using the four linearly independent solu- 
tions (ui,Vi),i=l-4, which satisfy the homogeneous system obtained from (2) 
by replacing the right hand side by zero [4,13]. Two of them are regular at 
r — > and the other two are regular at r — > oo. The latter are chosen with 
the asymptotic momenta k± = \/n ± E such that lmk± > which provides, 
in particular, the correct asymptotics for the particle scattering states at the 
IraE = axis. Such boundary conditions ensure the spatial localization of the 
nuclear ground state at ji < [3,4]. 

It is convenient to study the pairing effects in terms of the spectral distribu- 
tions over the eigenstates of the HFB Hamiltonian, i.e. in terms of the spatial 
integrals of the imaginary parts of the Gor'kov Green's functions. Here, for the 
given jl quantum numbers, we shall consider two partial spectral distributions: 
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The first and the second we shall call the v 2 - and w 2 -spectral density, re- 
spectively. The corresponding analytical expressions can be obtained by using 



the partial Green's functions which incorporate the generalized Wronskian of 
the four linearly independent solutions as given in Refs. [4,13]. We add that 
the expressions for the canonical-basis spectral distributions and the numeri- 
cal results with discretized continuum for a few tin isotopes can be found in 
Ref. [6]. 

In eq. (3), the spatial integral converges for R — * oo since it includes only 
the localized components oc v 2 . At E > |/x|, all HFB eigenstates are in the 
continuum and, since ImgY contains non-localized components oc u 2 , the in- 
tegration in eq. (4) is meaningful only up to a certain finite range R. One may 
present the two above expressions by one function dA^ FB /de if (3) and (4) 
are considered as functions of e = \x — E and e = \x + E, respectively, though 
both are for the same HFB eigenvalue E. We show in Figs. 1 and 2 the par- 
tial HFB neutron spectral densities for even and odd /, respectively, calculated 
for the drip-line nucleus 70 Ca, as functions of e with R=10 fm (lower panels) 
and also present the corresponding "HF" level densities (upper panels). The 
"HF" refers to the results obtained by solving eq. (2) but keeping only the 
diagonal part of Ti which emerges after solving the HFB problem (hji includes 
the contribution from pairing correlations). In other words, in these figures we 
show the level densities for the basis which diagonalizes the mean field part 
of the HFB Hamiltonian versus the spectral densities which arise when the 
pairing field A is present explicitly in the left matrix in eq. (2). We point out 
that the spherical volume of radius R=10 fm is chosen here as the physical 
volume in which the mean-field and pairing correlations are mostly active to 
produce the localized ground state, it is not a normalization box; within this 
volume we integrate the imaginary parts of the Green's functions which have 
the correct asymptotic behavior (the normal and anomalous densities, when 
solving the HFB problem, have been evaluated up to 25 fm, and, again, the 
wave functions are not zero at and beyond this radius, but the contribution 
from this far region to the bulk ground state properties is negligible). For the 
visualization of the 5-like peaks associated with the discrete "HF" states and 
very narrow resonances, a small imaginary part 7=1 keV has been added to 
the energy E in the Green's functions. The obtained chemical potential is 
small, \x = —149 keV; the spectral densities within the region 2yU < e < 
should vanish, and those which appear in Figs. 1 and 2 at these energies are 
due to the smearing parameter 7. 

Let us discuss first the peculiarities in the "HF" level densities. For a finite- 
depth potential well, such as our neutron HF potential, it can be shown that 
the partial single particle level density in the continuum within a spherical 
volume with radius R larger than the potential range is given by 
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Fig. 1. Neutron even-^ partial spectral densities in drip-line nucleus 70 Ca calculated 
by integrating the imaginary parts of the corresponding Green's functions up to 
R = 10 fm. A smearing parameter 7 = 1 keV is used to visualize the <5-like peaks 
and very narrow resonances. Upper panel: level densities for the "HF" potential; 
the e < part for the discrete spectrum with indicated single-particle levels, the 
e > part for continuum states. Lower panel: the HFB results; all states are in the 
continuum, the left part (e < 0) is for the localized v 2 -spectral density, the right 
(e > 0) for the non-normalizable u 2 -components (see text). 



where YJ t (kR) = (cos S ji)ji(kR) — (sin S ji)rii(kR), i — 1,1 ± 1; k — y/2me/h. 
This formula is valid for R big enough so that the radial wave function can be 
expressed through spherical Bessel and Neumann functions (ji and rii, respec- 
tively) and the phase shift Sji(e) can be calculated accurately. The resonance- 
like bumps seen in the upper panels in Figs. 1 and 2 at e > are associated with 
the term d5ji/de in (5). These bumps are imposed on a rather flat nonreso- 
nant continuum which weakly oscillates with an energy- dependent wavelength 
defined by A(kR) ~ 2ti. Without potential, when 8ji(e) = 0, eq. (5) gives the 
free-gas partial level density. Subtracting the latter from eq. (5), at kR ^> / 
one gets 
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Fig. 2. The same as in Fig. 1 but for the odd-£ states. 



Averaging this expression over small energy range Ae > nJ2h 2 e/mR 2 elimi- 
nates the oscillating cosine term, and one is left with the usual level density 
oc dSji/de associated with the potential well itself (see also Ref. [16]). How- 
ever, this procedure can not be applied at low energies, i.e. in the region 
near the continuum threshold which, in drip-line nuclei, is very close to the 
Fermi surface. This zone of the phase space is of importance since the ma- 
jor pairing effects are developed around the Fermi surface. At any finite R, 
when e — > 0, the leading term oc k 2l ~ x of the expansion of dSji/de in pow- 
ers of k cancels in eq. (5) due to the contribution coming from the term 
oc sm(2Sji)ji-i(kR)ni+i(kR). The resulting level density near the continuum 
threshold turns out to be proportional to Cjik 2l+1 , i.e. it has the free-gas be- 
havior, the coefficients Cji depend, however, on the potential well. For the 
s-wave level density, using the effective range approximation, we get from (5) 
and (6) an analytical expression: 



dNl 

2j 

de 



HF 



9 P 
— -R 3 

3n e 



4 



1-^ 
R 



1 _ lTcS 
2 R 



(7) 



where a is the scattering length and r c g is the effective range. In the curly 



brackets here is a factor which could give an enhancement of the HF level 
density as compared to the free-gas one. For our "HF" potential we found 
a = —14.8 fm and r e g- = 14.4 fm, which gives, with R = 10 fm, an enhancement 
factor of about 7. These parameters correspond to the virtual s-state at e = 
—51 keV which shows up as a sharp peak at very low positive e in the upper 
panel in Fig. 1. The width of this very asymmetric peak at half maximum is 
about 23 keV. At larger energies, one can see in this panel two partners of the 
o?-wave resonance which is split due to the spin-orbit potential, and also the 
#7/2 resonance whose gg/ 2 partner is in the discrete spectrum at e = —276 keV. 
Besides, a negative-parity h\\ii resonance can be seen in the upper panel in 
Fig. 2. The resonance characteristics and the energies e° of some bound "HF" 
states are listed in Table 1 (the first three columns). 

The graphical presentation of Fig. 1 and 2 has the advantage that the HFB 

Table 1. Characteristics" of some "HF" states and localized 
u 2 -components of the HFB resonances in 70 Ca 
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a Energies and widths are in MeV. 
b Virtual state (see text). 



spectral densities from the lower panels would be smoothly transformed to 
the corresponding upper "HF" ones if the A's in the HFB Hamiltonian of 
eq. (2) could be smoothly switched off keeping the "HF" part hji unchanged 
and defining \x by the particle number condition. In this way, one can easily 
trace the genetic origin of the irregularities seen in the HFB spectral densi- 
ties. In particular, the resonance contributions from the particle continuum 
to the localized v 2 -spectral density can be recognized. The characteristics of 
the picked-out localized v 2 -parts of the HFB resonances are listed in Table 1 
attributing them the genetic "HF" jl quantum numbers (the fourth and fifth 
columns). The widths of the "HF" resonance states and of the localized v 2 - 
components of the HFB resonances have been estimated from a 3-point no- 
background Breit-Wigner fit of the corresponding maxima. One can see that 
the discrete "HF" states get widths, i.e. they become the HFB resonances 
whose v 2 -components are shifted, as expected, to lower e. The HFB nonres- 
onant background can be also seen in the lower panels of Figs. 1 and 2. The 
non-localized u 2 - and localized v 2 -components of the HFB resonances are lo- 
cated symmetrically with respect to e = // (i.e., E — 0). We remark that all 
these belong to the same HFB eigenstate, however we do not discuss here the 
properties of the -u 2 -components due to lack of space. 

As seen in Fig. 1, the lgg/2 "HF" state at t\ = —276 keV becomes, with 
pairing, a HFB resonance whose localized v 2 -component lies well below the 
continuum threshold, at e ~ —A where A is the diagonal matrix element of 
the pairing field over the "HF" states near the Fermi surface. In our case, 
the energy gap 2A between the two components of this resonance is rather 
big, ~ 10 MeV. This could lead to the noticeable reduction of the tail of 
the neutron density distribution (the so-called pairing anti-halo effect [8,7]). 
However, in the drip-line nuclei, when |/x| <C A, near the threshold there 
appear strong localized v 2 -components from the resonant and non-resonant 
parts of the particle continuum with different jl. This is clearly seen in the 
lower panels of Figs. 1 and 2. 

The "HF" and HFB partial neutron occupation numbers Nji in 70 Ca are listed 
in Table 2. Also given in this table are the numbers of neutrons Nji(r > Ri) 
and Nji{r > R2) in the tail of p n beyond the nuclear surface; R± and R2 are 
defined by p n (Ri) = 0.01 fm~ 3 and p n (i?2) = 0.001 fm -3 , respectively. The 
choice of R\ roughly corresponds to the radius Rf + t/2 where t = 4(ln3)<i 
is the width of the diffuse zone approximated in the usual way by the Fermi 
function oc {1 + exp[(r — Rp)/^}' 1 . The calculated densities for Ca isotopes 
with A = 38 to 70 (a few of them are shown in Fig. 3) can be well described 
in the surface region by such a Fermi distribution. For 70 Ca, we found Rp = 
4.60 fm and d = 0.73 fm which gives R± = 6.2 fm (for comparison: in 48 Ca, 
Rf = 3.85 fm, d = 0.52 fm and R\ = 5.0 fm; thus, approaching the drip 
line, the size of the diffuse surface zone is significantly increasing). Using a 
larger radius R2, one can characterize the jl composition of the outmost tail 
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Fig. 3. Neutron and proton densities for some Ca isotopes resulting from the spher- 
ical HFB calculations with functional DF3 and gradient pairing force of eq.(l). The 
curves are marked with the mass number A. The two inserts are both for even- and 
odd-A, up to the drip-line nuclei 70 Ca and 55 Ca, respectively. Right upper insert: 
the number of neutrons in the "halo" region where the density does not exceed the 
indicated values (0.01, 0.003 and 0.001 fm -3 ). Left lower insert: the rms matter 
radii in comparison with "constant-pairing" HFB results (see text). 

of p n in the "halo" region. Table 2 shows that the pairing makes the Nji 
distribution more "uniform" than in the "HF" case, the major effect is due to 
the depletion of the discrete "HF" orbitals #9/2, /V/ 2 and / 5 / 2 and filling the 
resonance structures from the "HF" continuum (3si/ 2 , 2<i 5 / 2 , etc., see Table 1). 
The partial neutron numbers in the "halo" region are also distributed more 
uniformly, the contributions from the states with higher angular momenta, up 
to I — 8, are not negligible. Thus, in our calculations, the composition of the 
extended tail of the neutron density distribution in 70 Ca is determined by a 
collective pairing effect. 

The numbers of neutrons beyond the diffusive surface, iVhaio, calculated both 
for even and odd Ca isotopes with A = 38 to 70 are shown in the right insert 
in Fig. 3. Depending on the choice of the upper limit for the tail density, 
0.01, 0.003 or 0.001 fm" 3 , iV halo varies from 1.41 to 5.30, from 0.47 to 2.35 
and from 0.19 to 1.04, respectively; the relative increase of iVhaio for the latter 
choice reaches a factor of 5.5. Irregulaties seen in in the behavior of iVhaio are 
due to the odd-even effect. Odd isotopes are calculated with blocking [13], 
the heaviest one is predicted to be 55 Ca. The neutron density in this drip-line 
nucleus is also plotted in Fig. 3, the one- neutron "halo" tail of p Q in this case 
originates from the 2pi/ 2 blocked level with eA = —820 keV. 



The rms matter radii r m are shown in the left insert in Fig. 3. The staggering 
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Fig. 4. Binding energies E\, (lower curves, left scale) and two-neutron separation 
energies S^n (upper curves, right scale) in even Ca isotopes. Solid (dashed) lines 
connect the points resulting from the spherical HFB calculations with functional 
DF3 and gradient ("constant") pairing force (see text). Experimental data and 
those derived from systematics are from [17]. 

seen in the lighter Ca isotopes is generated by the gradient pairing force (1) 
whose parameters were chosen to reproduce both the neutron separation ener- 
gies S Q and anomalous behavior of charge radii [13]. One may suspect that this 
force yields too strong pairing near the drip line where the surface becomes 
more diffuse and the repulsive gradient term in (1) decreases; as the result 
we get A « 5 MeV for 70 Ca and A ~ 1.8 MeV around 44 Ca, the latter being 
in agreement with experiment. Unfortunately, the density dependence of the 
effective pairing force is not well known, and the extrapolation to the drip line 
with parametrization (1) may be questionable. Quite different results could be 
obtained with another kind of the pairing force. One can use, for example, in- 
stead of (1), a simple (5-interaction without p-dependence, with JF? = — O.66C0, 
which also reproduces well the experimental S n values, and calculate then the 
whole Ca chain. The matter radii obtained with such a "constant" force are 
shown by the dotted line in the left insert in Fig. 3, and the binding ener- 
gies together with two-neutron separation energies are presented in Fig. 4 in 
comparison with the results for the force (1). The behavior of r m versus A for 
"constant" T^ is relatively smooth, without noticeable staggering Fj the odd- 
and even- A drip-line nuclei occur to be 57 Ca and 70 Ca, respectively. Around 
these nuclei, weak irregularities can be seen in the behavior of r m . 70 Ca comes 
out with A = 0, i.e. pairing correlations do not appear for "constant" T^ 
in this nucleus; besides, the calculated S^ value is negative, —236 keV (see 



3 We remark that neither the observed evolution of charge radii in Ca nor the kink 
in Pb isotopes could be reproduced with this simple pairing force [13] which is thus 
rather unrealistic and serves here only as a demonstrative example. 
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Table 2. Partial neutron occupation numbers in 7U Ca 







Njl 


N jt (r 


> R 


i) a 


Nji(r 


>R2) b 


lj 


"HF' 


' HFB 


"HF" 




HFB 


"HF" 


HFB 








#1=6.10 fm 


Ri 


=6.15 fm 


#2=7.60 fm 


#2=8.15 fm 


Sl/2 


4 


4.740 


0.079 




0.501 


0.007 


0.176 


P3/2 


8 


7.861 


0.612 




0.711 


0.120 


0.146 


Pl/2 


4 


3.774 


0.408 




0.393 


0.104 


0.083 


d5/2 


6 


8.207 


0.094 




0.920 


0.006 


0.204 


d 3/2 


4 


4.725 


0.088 




0.439 


0.008 


0.103 


%2 


8 


8.012 


0.362 




0.554 


0.039 


0.077 


r 5/2 


6 


4.599 


0.589 




0.503 


0.129 


0.073 


§9/2 


10 


5.274 


1.267 




0.633 


0.294 


0.072 


§7/2 




0.877 






0.211 




0.035 


hll/2 




0.925 






0.191 




0.026 


n 9/2 




0.292 






0.097 




0.018 


il3/2 




0.321 






0.087 




0.014 


in/2 




0.136 






0.059 




0.011 


J15/2 




0.144 






0.052 




0.009 


J13/2 




0.054 






0.035 




0.008 


k 17/2 




0.040 






0.026 




0.006 


k 15/2 




0.019 






0.016 




0.005 


Xlj: 


50 


50.00 


3.500 




5.429 


0.707 


1.067 



a R\ is defined by the condition p n (R\) = 0.01 fm 3 
b #2 is denned by the condition p n (#2) = 0.001 fm -3 

Fig. 4), so 70 Ca could be a two-neutron drip-line double- magic emitter. As seen 
in Fig. 4, the binding energies E\> and ^n values are reproduced fairly well 
with both kinds of the pairing force, though a better description is obtained 
with the gradient force (1). Systematic deviations start at the mass numbers 
right beyond the measured region. Clearly, new experimental data and more 
extensive theoretical studies are needed to find still better parametrizations of 
the effective pairing force and improve the predictive power of self-consistent 
mean-field models. 

In conclusion: we have presented a numerically exact solution of the HFB 
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problem for the cutoff local energy- density functional. The Green's function 
formalism with physical boundary conditions has been used. The localization 
of the particle continuum in the drip line nucleus 70 Ca has been analyzed in 
terms of the partial spectral densities for the HFB eigenstates whose spectrum 
is fully continuous. The genetic relation to the "HF" (mean-field) partial level 
densities has been discussed, and the contributions from the resonant and 
nonresonant particle continuum to the localized ground state, as well as the 
shifting and broadening of the discrete hole orbitals have been clearly revealed. 
The composition of the extended neutron distribution in the "halo" region be- 
yond the diffuse surface has beeen analyzed in terms of the partial occupation 
numbers, and its formation has been shown to be a collective pairing effect. 
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